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General expressions of the electron yield in 57 Fe integral conversion electron Mossbauer spec- 
troscopy were derived depending on the angle of incidence of the 7-rays, on the source polarization 
and on the isotopic abundance of the source and the absober (blackness effects) using an expo- 
nential escape function of the electrons originating from all Mossbauer-resonance-related processes. 
The present approach provides a firm theoretical basis to determine the alignment and direction of 
^ ■ magnetization in the absorber. The intensity formulae were justified by least squares fits of a— 57 Fe 

^ ' spectral intensities measured in linearly and elliptically polarized source and absorber geometries. 

I The fits reproduce the experimentally set angles with high accuracy. Limits of the current approach 

and its relation to other, less complete treatments in the literature are discussed. 
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I. INTRODUCTION 



Mossbauer spectroscopy has been very successful in determination of the alignment and direction of magnetization 
of buried layers of Fe, the line intensities being dependent on the angle between the propagation direction of the 
X-ray beam from the source and the hyperfine field. Using linearly polarized source we may determine alignment, 
and, using circularly polarized radiation, also the sign of the magnetization component parallel to the propagation 
'-^J . direction^. For this reason polarized Mossbauer sources were prepared^. The so called filter technique^ provide 
almost single line circularly polarized radiation, resulting in relatively simple spectra, but uses only a small part of the 
O . source intensity, application of a split multi-line sourceii^, however, results in more complex spectra, but provides 
1^, ' higher intensity. 

Due to the relatively shallow escape depth of the electrons, conversion electron Mossbauer spectroscopy (CEMS) 
t-H ■ is used for suface and thin film investigations. Not even in the case of homogeneous magnetization are conventional 
CEMS experiments (in which the 7-rays are perpendicularly incident on the sample plane) are appropriate to determine 
. an arbitrary direction of magnetization. The relative intensities do not depend on the alignment of the magnetization 
within the sample plane, a frequent case with thin films, in which the shape anisotropy constraints the magnetization 
I , in the sample plane. 

The theoretical description of transmission Mossbauer spectroscopy is available for quite general cases. However, 
in a CEMS experiment, the electron transport from the Mossbauer nucleus to the sample surface is rather complex 
■ which renders the calculation of the electron yield difficult. There exist semi-empirical formulae (based on Monte- 
I Carlo simulations) providing the probability of the escape of an electron emitted with a given energy at a given 
depth below the surface^, which can be built into the numerical algorithms, but seemingly do not allow analytical 
solutions, which would provide more thorough insight. As a result, most spectroscopists, when evaluating CEM 
spectra, neglect the use of such escape functions and assume that the sample is very thin. In case of perpendicular 
incidence and non-enriched samples, this may be justified, unlike for slanted incidence and for absorbers enriched in 
the Mossbauer-active isotope, when polarization- and angle-dependent blackness effects are appreciable. Even from 
otherwise quite general treatments of multilayer structures, e.g£ discussion of polarization and blackness effects are 
missing. Papers that qualitatively elaborate such effects e.gi 7 -, leave room for improvement in regard of generality 
and mathematical thoroughness. The opposite extreme is the grazing incidence (below 1°) CEMS spectroscopy-^, for 
which the theory is well established, but, as we shall see, at such angles, a different approximation of the theory 
applies. Moreover, experimental feasibility of grazing incidence CEMS using radioactive sources is quite limited due 
to solid angle reasons, which enables the use of less than a fraction of w 10 of the source intensity. 

In the last decade several instruments and methods were devised in our laboratory for CEM polarimetry at arbitrary 
(except grazing) angle of incidence and linear as well as elliptical polarization ^ 10 i 1 1 . The analysis of the measured 
spectra with these setups lead us to recognize, that treatment of the blackness effect and polarization of the radiation 
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in the literature for conversion electron Mossbauer spectroscopy at low angle of incidence and of enriched samples in 
the resonant isotope is not satisfactory^. The present paper tries to fill this gap. Here we derive the CEM spectral 
intensity in case of polarized (Mossbauer) sources, when the 7-rays fall on a homogeneously magnetized sample at 
arbitrary (except for grazing) angle of incidence. We will show under what circumstances the escape function should 
be used and be neglected and the absorber consequently be considered thick and thin, respectively. 

II. CEM SPECTRA WITH THE POLARIZED 7- RAY INCIDENT ON THE ABSORBER AT AN 

ARBITRARY ANGLE 



In the following, concerning nuclear resonance we will follow the theoretical considerations and formalism used i: 
Here the case of a single, homogeneous layer will be considered. Generalization for multilayers is staightforward. 
The number of detected electrons created by a photon beam of energy E may be written as: 



,12 



D 



A/"ce(A e ) oc j Tr (p(z, E)N(z)a AI (z, £)) W(z)dz, (1) 


where D is the sample thickness, ~p{z) is the density matrix of the photons at a depth of z, N is the density of the 

scattering centres, a M is the part of the absorption cross section tensor belonging to Mossbauer effect, W(z) is the 
probability of a conversion electron be emitted following the absorption of an incoming photon at depth z and that it 
reaches the surface of the sample. In the case of differential CEMS (DCEMS) one should write W(z, E e ), where E e is 
the energy of the emitted electron, and W(z, E e ) gives the escape probability of an electron being emitted with energy 
Ej^. For brevity, unless it may lead to misunderstanding, we will omit the energy dependence in the formulae. 
The transfer operator T(z) defined by 

\iiz)) = i( Z )\m) (2) 

determines the wave at depth z from the wave \ip(z = 0)) at the surface. Accordingly, for the density matrix we get 

^(z)=T(z)^(0)f + (z), (3) 
where T denotes the hermitian adjoint operator. After some permutation and substitution we may write 



Nce(D, E) oc Tr (p(0, E)J(D, E)j , where 



D 



J(D,E)= J T + (z,E)a M {z,E)T(z,E)N(z)W(z)dz. 



(4) 





For a single homogeneous layer T(z) has the following forrn^: 

T(z) = sinh fifcomz sin $) m" (T — R) + cosh {ik^mz sin #) (T + R), (5) 

where fco is the wave number in vacuum, 1? is the glancing angle of the incoming X-ray beam, I denotes the unity 

matrix, | '(/'Reflected (0)) = R |^(0)) is the reflected wave (with R being the reflectivity tensor) at the surface of the 
sample and we introduced the tenzor quantity TTT with the definition 



m= \A+-^X, (6) 
V sin v 

where x is the dielectric susceptibility tensor. For multilayers the resulting transfer matrix is the product of single 
layer T matrices of form (JSJ). For large angles of incidence, ■& the reflected beam for X-rays is weak, i.e. R is small as 
compared to the unity matrix. Consequently 

T(z) = sinh (ikomz sin j?) m 1 + cosh (ik^mz sin #) . (7) 
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For nearly perpendicular incidence ($ ~ ^) the matrix m will be the index of refraction matrix n = y I + \ ~ I + \x- 
Since the dielectric susceptibility for X-rays is always small, the inverse of the index of refraction matrix is also close 
to the unity matrix, n" ps I , which leads to the approximation given in^, namely: 



T(z) = e ifconz « e ikoZ e^ koXZ . (8) 
For small angles of incidence approximation ((8|) is not valid. We may expand, however, m and m 1 into Taylor-series 

= t 1 = 1 =2 , 

;X - n . 4 a x +■■■ 



2sin 2 i? 8sin 4 tf 
_L - 3_ 

2sin 2 tf X 8sin 4 tf 



m T ;— • — r^x • ••• (9) 



In first order we get 

T(z) w c i ' £o2sin ''e i ^ of ii^ + sinh(ifc rrTzsintf) (rrT 1 -T) . (10) 



Higher order terms are negligible if the eigenvalues of \ fulfill the condition \\%\ < sin 2 For 57 Fe on resonance 
(IXil « 0.01), therefore the first order approximation is valid down to about & w 4 — 5°. Unless for very small angles, 

m « T and therefore TrT 1 ps T. Since in the following we will not consider grazing angles, we will neglect the second 
term sinh(. . . )(m — T) on the right side of (jTUJ) . 

The absorption cross-section tensor for a single atom (molecule, etc.) can be defined as&i£ 

= « - J^T -X= + (*)) ■ (11) 



2V(z) V w V V 2iV(z) 

Near the nuclear resonance, the electronic contribution to the susceptibility is negligible compared to the nuclear one, 
therefore we may write: 

Substituting (Ti"2"|) and the first term of pop into (J3J), the CEMS intensity matrix J for a homogeneous layer becomes: 

J(Z>) = - J e-^^ii^^ (f e i ^«SI^PF(z)dz 



= -sim?y VF(z)— ^e"^ feo ^sn3 e'^^snjj dz. (13) 

Introducing the s = iifcoz/ sini? notation and using the Baker - Campbell - Hausdorff formula^ the product of the two 
matrix exponentials can be written as 

e _s ^ + e s ^ = e s (^"^ + W s2 ^ + 'fl + A s3 [^-^ + .[^ + <!]]+■■■. (14) 

For the Ml nuclear transitions, like in 57 Fe, for almost all practical cases the x + ] = condition is fulfilled with a 
good approximation-ii. Otherwise one has to estimate the contribution of higher order commutator terms in (I14p . In 

the following we will assume that [^,x + ] = is fulfilled, therefore we have 



i{D,E) =-wa.<ej W ^§; (e- NaM(E) ^^) dz. 



(15) 



Not regarding the number of approximations we had to make in order to obtain cq. (|15|) . the result may seem to be 
commonsensical: the sm r & term in front of the integral simply takes the reduced cross section of a slanted sample into 
account. The number of resonantly absorbed photons at depth z is given by the derivative term inside the integral. 
For small angles the photons travel a longer path z / sin $ to reach depth z and therefore they get absorbed more likely 
near the surface as compared to the case of perpendicular incidence. However, one should keep in mind, that eq. (JTS 
is only valid for not too small angles, when the linear approximation in eq. (JTUJ) is valid. 
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Until now we did not take an essential feature of the Mossbauer effect into account, namely the moving source. The 
density matrix of the 7-ray beam emitted by the Mossbauer source that moves with velocity v can be written as 

where p^" 31 characterizes the intensity ratio and polarization of the i-th line at position Ei in the energy spectrum, 
Ci{E 1 v) is the (Lorentzian) shape function of the i-th line, T is the line width and c is the speed of light. As a function 
of the Mossbauer drive velocity the measured CEMS intensity is 



/cems(A«) = isI^Tr (Pi°J d(E,v)J(D,E)dEj , (17) 

where fs is Mossbauer-Lamb factor of the source and I is the total intensity of the source incident on the sample. 

In order to calculate the CEMS spectra the knowledge of the electron escape function W(z) is necessary. Liljcquist 
derived semi-empirical formulae based on his Monte-Carlo simulations^^. Instead of depth, he uses the equivalent 
depth in iron defined as 

t = -?-z, (18) 

QFc 

where g is the density of the sample and ppe the density of natural iron at ambient conditions. For a homogeneous 
layer, the escape function can be well approximated with an exponential, 

W(z) = Woe'^ 1 = W c~ Z ~^ = W c~C(e) ( (19) 

where Wq expressing the probability of emission of a conversion electron after photon absorption at an iron-equivalent 
depth of r w 540 FeA to reach the surface of the sample. The more general escape function are derived in appendix 
lA"]) . The exponential approximation is advantageous in the parameter fitting routines, since the code based on this 
is by two orders of magnitude faster than the one based on exact formulae, while providing a good approximation. 
Moreover, the exponential approximation is easier to handle analytically in some special cases. Indeed, with (|15[) we 
have 



/cEMs(A«) = a]T> Th° l J C i (E,v)CNW M (E)smd 



(NW M (E) +Tsini?) (T-e ° 



dE 



(20) 



where a = WofsX. 

First let us consider the unpolarized case. Applying the formula of effective thickness as t — g/ gp c D = N/N-p e D, 
and similarly r = A/Ap e C an d the absorption coefficient (for natural iron at ambient conditions) as a = N-p c a, the 
exponent in eq. (|20p can be transformed as follows: 



D 



( Na(E) 1 
It may be convenient to interpret the quantity 



-ext- 



ra. + sin i9 



1. sint? 



(21) 



t'(E,ti) = t 



ra(E) + sintf 
ra(E) sin?9 



(22) 



as a corrected effective thickness^. 

In the polarised case for small angles (sini? 0) we have (Na (E) I (Na (E) + I sin^ J « I and therefore 



Icems(D,v) = asin-d 



=l E Tr (^ 0l / k {E, v) c"^ M ^ & dE 



(23) 
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which formally corresponds to the thick sample case in absorption Mossbauer spectroscopy, as it was conjectured ini. 
We should not forget, however, that beside the above mentioned constraints in deriving (Ti~S"|) and (|2"0|) , this approxi- 
mation is not valid for grazing angles, when the reflections due to the electronic contribution to the susceptibility are 
not negligible. 

At nearly perpendicular incidence, when sini? m 1 we expand the exponent in Taylor series up to second order: 

ft - e -»0v*"<*>+iD) „ D (nW M (E) + If) + (Naf(E) + £f + ... (24) 
In case the second term is negligible compared to the first one, i.e. 

l -D 2 (Na™(E) + ^ « D (Naf(E) + , (25) 
where Oi is the z-th eigenvalue of cf, the intensity expression 



JcemsCAtO = a^D^Tr (rf* J Na M (E)Ci(E,v)dEj , 



(26) 



formally corresponds to the thin sample case in Mossbauer spectroscopy, as discussed by£. 
From the above, the condition of the thin absorber limit in perpendicular incidence is 

^aT^TT' (27) 

where a^ ax is the maximal eigenvalue of the absorption cross section tensor. 

This condition means, that D is the characteristic information depth of the surface of the sample for CEM 
spectroscopy; no significant number of electrons will reach the sample surface from deeper layers of the sample and/or 
the photons will not reach a larger perpendicular depth in the sample. 

For multilayers we may follow the same steps, the main difference and question will be the electron escape function, 
which may become quite complex. The conversion electrons emitted from atoms in different layers should be summed 
up, taking into account the proper escape function(s), which probably may be determined as some good examples of 6 
show, but this is out of the scope of this paper. 



III. THICK SOURCE 



Blackness effects in Mossbauer spectroscopy are not limited to the absorber. Due to the possible self-absorption 
of the emitted photons, additional broadening may arisei^ or, in case of a split multi-line source, the relative line 
intensities of the emitted lines of the Mossbauer source may considerably be modified, provided that 7-emission occurs 
relatively deep below the surface of the source and the substrate of the source contains large amount of the resonant 
isotope. In the commonly used single line sources these effects are avoided by letting the resonant isotope diffuse only 
into a shallow surface layer of the substrate. When using such sources in perpendicular emission geometry, blackness 
effect rarely occur. However, the relative intensities of a split multiline source of a-Fe substrate at slanted exit of the 
7-rays, which is studied experimentally and discussed below arc considerably modified by the blackness effects. In 
the following we discuss the case of the thick source. 

Following the same logic as in deriving ([3]) the density matrix of the photons emitted by the Mossbauer source may 
be written as 

p(E) oc J 1 + (D S - z, E)W M (E)1(D S - z, E)Af s (z)dz, (28) 


where D s is the physical thickness of the substrate of the source and A/" s (z) is the number of the nuclei emitting the 
resonant photons at depth z. The Mossbauer sources are commonly prepared by introducing the radioactive atoms 
into the host material by diffusion. Therefore it is correct to assume that Af s (z) follows a diffusion profile, i.e. an 
expression containing an error function. 
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FIG. 1: Results of linear polarimetric measurements (circles) and their fits (lines) with parallel and perpendicular source and 
sample magnetizations. The X-ray beam is perpendicular on the plane of the source and the sample. The magnetizations are 
in plane. 

7500 
6000 
4500 



9600 
8400 
7200 
6000 

-15 -10 -5 5 10 15 

velocity (mm-s 1 ) 

FIG. 2: Results of circular polarimetric measurements (circles) and their fits (lines) with parallel and antiparallel source and 
sample magnetizations. The X-ray beam is enclosing § — 10° with the plane of the source and i? = 5° with the sample. The 
external magnetic fields are parallel to the X-ray beam. 

IV. EXPERIMENTS AND THEORY 

CEM polarimetric measurements were performed using a-iron foils, a natural Fe foil of 15 /urn thickness as the 
substrate of the 57 Co(a-Fe) source and a 57 Fe foil of 20 fxm thickness as the sample. A detailed description of these 
experiments was published inia. The experimental spectra are the same here as ir>», but the spectra are evaluated 
using a different model function. The emphasis here is being on the verification of the theory derived in the present 
paper. 

Linear polarimetric measurements were performed in perpendicular incidence geometry, both the source and sample 
was magnetized perpendicular to the propagation direction. The results of such two measurements are shown in Fig. Ml 
where the magnetic fields on the source and sample are parallel and perpendicular to each other. 

Elliptical (close to circular) polarimetric measurements were performed with a slanted source and a slanted sample, 
both magnetized in their planes the magnetization vectors and the wave vector laying in the same plane. The sample 
planes enclosed d = 10° (source) and ■& — 5° (sample) with the propagation direction. The results of such two CEM 
experiments are shown in Fig. [5] where the external magnetic fields (230 mT) on the source and sample are parallel 
and antiparallel to each other. 

Having implemented the theory outlined in the former sections, measured spectra were fitted using the general fitting 
program FitSuite developed earlier—. The positions and intensities of the Mossbauer lines (i.e. the contribution of 
the resonant absorption to <t(£7)) were calculated solving the eigenvalue problem belonging to the corresponding 
hypcrfinc Hamiltonian. The theory descibed and implemented in EFFIsa were adopted and modified to include 
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the theoretical considerations presented above, namely eq. (|20]) . 

TABLE I: Essential fitted parameters, linear polarimetry. $b,... are the angles enclosed by the hyperfine field and the z axis 
(propagation direction of the X-rays). v?b,... are the angles are the azumithal angles measured around the z axis from the x 
axis being in plane of incidence. 



Parameter 


Fitted 


Expected 


|B| 


33.04 T ± 0.02 T 




$ B, source 


86.0° ± 0.9° 


90° 


V-'B, source 


0° (fixed) 




Parallel case: 






^B, sample 


108.9° ± 0.7° 


90° 


V-'B, sample 


4.2° ± 1.7° 


0° 


Perpendicular case: 






^B, sample 


90.1° ± 4.1° 


90° 


V-'B, sample 


91.2° ± 4.8° 


90° 



In linear polarimetric case the X-ray beam is nearly perpendicular to the source and sample planes, therefore no 
thickness and angle dependence is expected according to (|2"6"j) . The magnitude and orientations of the hyperfine fields 
(see TablelU) were found to agree quite well with the experimental values, set by the direction of the external magnetic 
fields. The fits and their x 2 statistics are acceptable. In agreement with this according to our simulations change of 
the thickness give rise to intensity change, but the shape of the spectra remains unchanged. 

TABLE II: Essential fitted parameters, circular polarimetry. #b,... and <Pb,. . . are defined in the same way as in Table |T] 



Parameter 


Fitted 


Expected 


Source slant angle 


6.53° ± 0.03° 


10° 


Sample slant angle 


4.936° ± 0.003° 


5° 


|B| 


32.9667 T ± 0.0004 T 




^B, source 


8.2° ± 0.6° 


10° 


V-'B, source 


0.3° ± 0.08° 


o° 


Parallel case: 






^B, sample 


13.8° ± 0.3° 


5° 


V^B, sample 


0.16° ± 0.09° 


0° 


Antiparallel case: 






^B, sample 


193.9° ± 0.2° 


185° 


V^B, sample 


0.2° ± 0.36° 


0° 



This is not the case for circular polarimetric experiments, where the blackness effect has an important role and 
therefore the intensity ratios strongly depend on the thickness and slant angle, as it can be seen in Figs. [5][31 Without 
this effect, the intensity ratios of the spectra of Fig. [2] could not be explained either. Although the features of the 
measured spectra are reproduced relatively well in this case, in the antiparallel case the figure of merit of the fit 
worse, and the fitted hyperfine field orientations (Table [TTJ) slightly differ from the actual experimental values. These 
discrepancies are probably due to fact, that the reorientation by the small external field was not completely succcsful. 

In fitting the circular polarimetric spectra at an angle of incidence of t? = 5°, when the perpendicular information 
depth D is much lower, we had to take into account an additional singlet site with isomer shift (0.11 ± 0.02) rrrm-s -1 , 
effective thickness (3000 ± 100) Fe 57 A with Lorcntzian broadening of (5.8 ± 0.6) T nat (r nat =0.097 mm-s -1 is the width 
of the natural line broadening) which is most probably due to the iron oxide on the surface of the absorber. In case 
of linear polarimetric measurements this component did not appear, since in pependicular incidence its contribution 
to the measured spectrum was much weaker. 

V. SUMMARY 

We have found that the CEM spectra may be regarded, as Mossbauer spectra: 
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FIG. 3: Demonstration of blackness effect because of thickness dependence. Results of circular polarimetric simulations with 
antiparallel source and sample magnetizations for different sample thicknesses. Other parameters are the same as in the 
corresponding fitted spectrum plotted in Fig. [2] 



• of a thick sample, if the X-rays coming from the source fall on the sample under small angle, which is still 
greater than what we usually call grazing incidence; 

• of a thin sample, if the X-rays coming from the source fall on the sample perpendicularly and the condition Q27|) 
is fulfilled. 

These statements were verified for electron escape functions of a single exponential function, in cases where this 
approximation is inappropriate these statements may not be valid either. Furthermore, we may not regard the CEM 
spectra neither as a thin nor a thick Mossbauer spectrum for angles between the two limiting cases. 

We fitted the experimental results according to the model based on considerations of this paper using the program 
FitSuite. 



APPENDIX A: ELECTRON ESCAPE FUNCTION 



The electrons are emitted because of several processes after irradiation by X-ray, which all should be taken into 
account in the escape functions. According to^ the electron escape probability from depth z can be written as: 



W(z) =^2ciWi{z), (Al) 



where Cj gives the contribution of the i-th process, for which the escape function is Wi(z). The escape function, using 
as argument the equivalent iron depth t instead of z (definition is given in (fT5|) ) for K, L-M and Auger electrons, has 
the form: 

w ( _ t \ _ { 0.74 - 2.7x + 2.5x 2 if s$ x < 0.55 (M) 
1 \ ~ rf ethe ) ~ I otherwise ' 
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FIG. 4: Demonstration of blackness effect because of angular dependence. Results of circular polarimetric simulations with 
antiparallel source and sample magnetizations for 40 nm thick samples slanted in different angles. Other parameters are the 
same as in the corresponding fitted spectrum plotted in Fig. [2] 



where rf ethe is the Bethe range (which is on average the total distance, which an electron travels until it is stopped 
down). Alternative formulae for K, L-M and Auger electrons based on Gaussian functions are available i n 18 ' 23 : 



W(t,E e ) = A(Z,E e )e 



9R(Z,E C 



(A3) 



where Z is the atomic number and E e the energy with which the electron was emitted from the atom 
details regarding the functions A(Z,E e ) and R(Z,E e ) seei 8 -. To get the alternative Wi{t) functions of 
to integrate W(t, E e ) according to the electron energies in the corresponding ranges too. Therefore 
appropriate for DCEMS problems, for which (|A2[) is completely inappropriate. 
For photo-electrons the escape function has the form: 



For further 
2]) we have 
3l) is more 



Wi(t) 



„cff 



1-0.5[E 2 GM) 
0.5 [E a (/ii(t- 



rf))-E 2 ( Mi i)] 



OsCfsC rf 
rf < t 



(A4) 



where ^{x) is the second order exponential integral function defined by E2(.t) = h 2 e xh h = J Q e " 
the so called effective range, fii is absorption coefficient for the conversion photons, Uooi is a normalization constant. 
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